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Merging the functionalities of computation process- 
ing and data storage in one single material plays a 
crucial role in the burgeoning field of spintronics[l, 2J. 
One of the promising candidates is the diluted mag- 
netic semiconductor 

Hi] 

such as (Ga.Mn)As, where 
itinerant carriers in the semiconducting bands mediate 
the effective exchange coupling between local mag- 
netic moments of the transition metals and lead to 
ferromagnetic ground state. Another approach, now 
growing into a subfield under the name of spin Hall 
effect, proposes to generate spin currents by electric 
gates via the spin-orbital interactions and has been 
verified in experiments. 

Here we would like to propose an alternative non- 
magnetic resolution, making use of the subtle interplay 
between Coulomb interactions and lattice topology at 
nanoscale. By both analytic weak-coupling and nu- 
merical density matrix renormalization-group methods, 
we show that novel itinerant ferromagnetism appears, 
upon appropriate doping, in (hydrogenated) armchair 
nanographite ribbons and carbon nanotubes. Inspired 
by recent experimental breakthrough in single-layer 
graphite[5, 6J, the metallic and ferromagnetic state 
we found here may have potential applications in spin- 
tronics at nanoscale. 

In this Letter, we study the ground state properties of 
the doped armchair nanographite ribbons (NGRs) and 
carbon nanotubes (CNTs) with short-range interactions. 
FigureHshows the lattice structure of the armchair NGR 
with open edges. Even though only tt electrons are ac- 
tive in low-energy, the honeycomb structure at nanoscale 
gives rise to both localized Wannier and itinerant Bloch 
orbitals. Furthermore, we show that these localized or- 
bitals form flat bands with zero velocity. Electron corre- 
lations in the flat band generate intrinsic magnetic mo- 
ments (see below) and the itinerant Bloch electrons medi- 
ate ferromagnetic exchange coupling among them. Thus, 
the magnetic and transport properties can be engineered 
by changing the width of the NGR and the electric field, 
opening up the possibility for future spintronics applica- 
tions at the nanoscale [ll Isl Wi IToj . 




FIG. 1: Armchair NGR of L x — 4 and L y = 5. Open edges 
are present at y = 1 and L y . The solid circles and squares 
represent sublattice A and B, respectively. The open circles 
show the amplitudes of the Wannier orbital of E = t at x — 2, 
with opposite signs indicated by light /dark shades. 



Besides its possible applications to nanospintronics, 
the ferromagnetic ground state in the doped armchair 
NGR also reveals a new type of flat-band ferromag- 
netisrm different from those proposed by Mielke and 
Tasaki pnLll^ . The key to Mielke- Tasaki ferromagnetism 
is the presence of local Wannier orbitals in the flat band. 
If the so-called local connectivity condition for adjacent 
Wannier orbitals is achieved, the finite overlaps gener- 
ate exchange coupling among these orbitals, leading to 
flat-band ferromagnetism. However, the Wannier orbital 
(shown in Fig. ^| in armchair NGR has zero overlap with 
its adjacent neighbors. As a consequence, the flat band 
alone only accounts for the existence of the magnetic mo- 
ments. As will become clear later, the ferromagnetic or- 
der sets in only when additional gapless itinerant carriers 
are present. 

To understand this new type of flat-band ferromag- 
netism in the armchair NGR, we start with the Hubbard 
Hamiltonian, H = H t + Hu, 



H = -t 



E 

(r,r'),< 



[c+(r)c a (r') + H.c] 
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where t is the hopping amplitude on the honeycomb net- 
work, U > is the on-site repulsion, a =|, J, is the spin 
index, r = (x,y), and Y]/ r r i\ is taken only for nearest- 
neighbor (NN) bonds. Non-uniform hopping due to lat- 
tice distortions and short-range interactions with general 
profiles will be discussed at the end. The values of t 
reported in the literature^! Q, range from 2.4-2.7 



2 



eV for CNTs, while t ~ 3 eV is typical in graphites. 
Although an accurate value of U is not yet known in 
nanographite systems, the estimate from polyacetylene, 
U ~ 6-10 eVlialla, might serve as a reasonable guess. 
Thus, we expect that U/t ~ 0(1) in nanographite sys- 
tems. However, to gain physical insight, it is helpful to 
study the weak coupling limit first. 

To simplify the problem, let us take L x — > oo here. 
In weak coupling, U/t <C 1, it is natural to concentrate 
on the hopping Hamiltonian H t first. Its band structure, 
shown in Fig. |2 can be obtained analytically by partial 
Fourier transformation and the supersymmetric (SUSY) 
algebra 18, 19]. Note that the flat bands located at E = 
±t appear only for odd L y . The local Wannier orbital at 
x = xq in the flat band takes the simple form, $ Fa:o (r) = 



(a) ^=5 



(b) L v =6 



$x,x ®F(y), where 



Ml/) 



sin(7n//2) 
=Fsin(7ry/2) 



(2) 



and the subscripts A/B label the two sublattices in 
Fig. ^ The =F signs stand for the flat bands at E = ±t. 
Repeatedly applying the lattice displacement operator T x 
on one Wannier orbital, all orbitals at different locations 
can be constructed. Since [T x ,H t ] = 0, all the orbitals 
share the same energy and form a flat band. 

At first sight, it is rather counterintuitive that the local 
Wannier orbital cannot move around by quantum hop- 
ping. The static nature is due to perfect destructive 
quantum interferences which make hopping amplitudes 
from different sites cancel each other. Furthermore, since 
the wave function is not zero only at odd y coordinates 
(see Fig. [TJ, it is clear that the adjacent orbitals have 
zero overlap. Thus, the Mielke-Tasaki mechanism does 
not work to couple neighboring orbitals magnetically. 

In the following, we concentrate on the flat-band 
regime E — t for the NGR with odd L y . To account for 
the low-energy physics correctly, it is necessary to include 
the dispersive bands intersecting the flat band at Fermi 
points k m = 2m7r/ (L y + 1), with m — 1,2,..., \{L y — 1). 
Applying the SUSY technique again, the Bloch wave 
functions of the itinerant carriers at the Fermi points are 



obtained as $ Pm (r) 



c $ Pm (y), where 



1 



e iPk m S A ( v ) s i n (fc m 7//2) 



(3) 



5 A/B (y) are the phase factors resulting from the honey- 
comb structure, and P = R/ L = ± denotes the chiral- 
ity. From the Fermi momenta k m , the lower and up- 
per bounds of the doping rate x = (n) — 1 for the flat- 
band regime are obtained as x m - m — ^2 m k rn / '(j{L y ) = 



1 



) and x Fl 



-, respectively. 



The effective low-energy theory can be derived from 
the chiral decomposition of the lattice fermion operator, 



c a (r) ~ <$> F {y)i) Fa {x) + 



°$ Pm {y)^ Pma {x), (4) 




k„/7t 1 



FIG. 2: Band structure for the infinite armchair NGRs of 
different width L y = 5, 6, 7, 8. 



where ip Fa (x) and ip Pma (x) are the fermion- field opera- 
tors for the flat-band and chiral states respectively. Since 
the density of states is divergent for the flat band, the 
itinerant chiral fields can be dropped to the leading or- 
der approximation . The effective Hamiltonian, keeping 
the flat band only, is rather simple, 



U E ^2 n Fl( X ) n Fl( X ), 



(5) 



where U E = U/(L y + 1) and n Fa (x) — i/j Fa (x)ip Fa (x). 
Since the kinetic energy is quenched in the flat band, the 
ground state contains no quantum fluctuations. To avoid 
the cost of U E , the ground state consists of the maximum 
number of singly-occupied Wannier orbitals, which leads 
to local magnetic moments. Since there is no overlap 
between adjacent orbitals, these magnetic moments are 
free and give rise to a large ground-state degeneracy. 

To lift the large degeneracy of the ground state, in- 
teraction with the itinerant carriers must be included. 
The effective Hamiltonian can be derived from the com- 
plete chiral-field decomposition in Eq. The dominant 
terms in the spin sector are the exchange interactions 
which couple the local moments and the itinerant carri- 
ers, 



H j — S J m S f (ie) • S m (ir), 



(6) 



where S F (x) = ^ip Fa (x)(T aa 'ip Fa >(x) is the spin density 
operator for the local moment (<r QQ < is the Pauli matrix) 
and S m (x) = ^ p S Pm (x) = £) P \il> Pma (x)cr aal ip Pma > (x) 
for itinerant carriers in each band. The coupling strength 
is computed from the exchange integral, 



J — 2 



$Ur>)$ Pm (r')V(r> - r" )<&*>")<!>,>") 



a y=odd " 

which is ferromagnetic for repulsive interaction [/ > 0. 

The exchange coupling J m tends to align the local mo- 
ments from the flat band because it does not cause any 
extra kinetic energy. The ferromagnctically aligned mo- 
ments act back on the itinerant carriers and induce fi- 
nite polarization in the itinerant spins. The interacting 
Hamiltonian H E + H , therefore shows interesting two- 
step flat-band ferromagnetism - the flat band gives rise 
to local moments without direct exchange coupling, while 
the presence of gapless itinerant carriers mediates the fer- 
romagnetic order. The significant feature of the armchair 
NGR/CNT is that, even within the one-orbital Hubbard 
model without any magnetic impurity nor additional lo- 
calized levels, the electronic correlations give rise to both 
the local moments and itinerant carriers simultaneously 
due to the peculiar topology of the honeycomb lattice. 

Now we turn to armchair NGRs with finite L x . If one 
imposes the periodic boundary conditions along the x- 
axis, the system becomes a short segment of armchair 
CNT. In this case, only when the quantized momenta 
k x = 2mr/L x coincide with the Fermi points k m — 
±2nm/ (L y + 1), the gapless itinerant carriers are present 
and the ferromagnetic ground state is realized. For 
open boundary conditions with finite L Xl the situation 
is much more complicated; k x is no longer good quan- 
tum number and the band structure can be deformed 
by finite L x . However, from numerical calculations for 
the tight-binding Hamiltonian H t , we have found that 
the energy spectrum around the flat-band level E = ±t 
is not affected by finite L x and can be accurately ap- 
proximated by the quantization rule k x = ln/(L x + 1) 
(I = 1,2, ...,L X ). When the quantized momenta k x co- 
incide with the Fermi points k m , ferromagnetism sets in 
with the help of these gapless itinerant carriers. There- 
fore, depending on specific choice of L x , the ground state 
of the armchair NGR/CNT in the flat-band regime can 
be ferromagnetic or Curie-like paramagnetic. 

To check whether the weak-coupling picture survives 
for the realistic coupling regime U/t ~ 0(1), we em- 
ploy the density-matrix renormalization group (DMRG) 
method[20l| to study the system in intermediate-coupling 
regime U/t = 4. To look for higher-spin ground state, we 
perform the calculation in a subspace preserving total 
spin of the system by using the SU(2)-symmetric density 
matrix. We compute the energy differences between the 
higher-spin ground states and lowest-spin ones defined 
as, 

A(N e , SS) = E Q (N e ,S + SS) - E (N e , S ), (8) 

where N e denoting the number of doped electrons from 
half filling (n) = 1 and E n (N e ,S) is the lowest energy 
in the subspace characterized by N e and total spin S. 
The lowest spin So is defined as So = for even N e and 
So = 1/2 for odd N e . The calculation is performed for 
the system with L y = 5 and L x = 2, 3. The number of 
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(a) L x =2, L y =5, U/t=4 (b) L x =3, L y =5, U/t=4 




FIG. 3: Doping-rate dependence of the energy differences be- 
tween higher-spin and lowest-spin states in the L y — 5 arm- 
chair NGRs with (a) L x = 2 and (b) L x — 3. The squares, 
triangles, and diamonds represent the data for SS — 1, 2, and 
3, respectively. Inset in (a): Total spin of the ground state, 
S S s, for L x — 2 and L y — 5 armchair NGR at different doping. 



(a) (b) 




FIG. 4: Ground-state distribution of (a) charge density 
(n(r)) — 1 and (b) spin polarization (S z (r)) for L x = 2 and 
L y — 5 armchair NGR with (N e = 5, SS = 1). The values of 
(n(r)} — 1 and (S iz (r)} are positive at every site, and the areas 
of the open circles are proportional to their amplitudes. 



kept states is up to 1200, typically corresponding to 400- 
500 SU(2) multiplets. The truncation error is typically 
10~ 5 for L x = 2 and 10~ 4 for L x = 3, and the results are 
extrapolated to the limit of zero truncation error. 

For L x = 2, the quantized momenta k x = In/ 3 co- 
incide with the Fermi points of the chiral fields, and 
therefore, the weak-coupling theory predicts the ferro- 
magnetic ground state. As shown in Fig. [21 we do 
find numerically the ferromagnetic (higher-spin) ground 
state in the flat-band regime (N e — 5,6,7). The nu- 
merical values of A for these higher-spin ground states 
are A(5,<55 = 1) = -0.036*, A{6, 5S = 2) = -0.074* 
and A(7, SS = 1) = — 0.035i, respectively. In addition, 
the local charge density (n(r)) — 1 and spin polarization 
(S z (r)) in the ground state are also calculated, shown in 
Fig. While the charge distribution is almost uniform, 
it is rather remarkable that the spin density has a simi- 
lar profile to the Wannier orbitals in the weak-coupling 
limit. This indicates that the physical picture developed 
in weak coupling still works rather well even for U/t = 4. 

For L x = 3, on the other hand, the approximative 
quantization rule gives k x = In /A, so that no gapless 
itinerant carriers in weak-coupling theory are available. 
In the numerical calculation, we find that in the flat-band 
regime, N e = 7,8,9,10, and 11, the higher-spin ground 
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state only shows up for iV e = 8, 10. The numerical val- 
ues are A(8, SS = 2) = -0.006*, A(10, SS = 1) = -0.02< 
- relatively smaller than those in L x = 2 case. Since 
the energy differences arc very small, it is quite natural 
to consider that the lowest- and higher-spin states are al- 
most degenerate, suggesting a paramagnetic ground state 
with almost free moments. We therefore believe that the 
numerical results support the novel finite-size effect pre- 
dicted by the weak-coupling analysis. 

Before closing, it is interesting to ask how robust the 
flat-band ferromagnetism is beyond the Hubbard Hamil- 
tonian. Note that the flat band in armchair NGR is 
unique in several ways. First of all, to achieve flat bands 
in other systems, one needs to fine-tune several parame- 
ters while it is not necessary for armchair NGR - when 
the transverse width L y is odd, the presence of the flat 
bands is guaranteed by the topology of the lattice struc- 
ture. In addition, in practical graphite network, the hop- 
ping along horizontal and titled bonds is expected to be 
slightly different. It is straightforward to show that this 
deviation only shifts the flat band from E = ±t but the 
bands remain flat. The profile of the short-range inter- 
action does not seem to do much harm either. In the 
presence of the NN interaction V, the exchange coupling 
in Eq. © becomes 

J m = —^—r(U -V cos k m ). (9) 

(Ly + I) 

Since V < U is often expected, the exchange coupling is 
still ferromagnetic and the picture does not change. 
Recently NGRs have been successfully fabricated pi) 



and may bring up the possibility to realize this new kind 
of flat-band ferromagnetism. Meanwhile, the develop- 
ing nanotechnology opens up another window to probe 
the interesting phenomena in a nanojunction by crossing 
a semiconducting CNT on top of a single-wall armchair 
one. Applying gate voltages in the semiconducting nan- 
otube to deplete electrons locally, it is possible to produce 
the desired spatial profile of the density so that some 
regimes are in the flat-band ferromagnetic/paramagnetic 
phase. The key issue is whether it is feasible to achieve 
sharp density profile by gate voltage alone. The most 
promising approach is provided by recent breakthrough 
in the fabrication of single-layer graphite 5, 6] with high 
mobility. Through the four-probe transport measure- 
ment, it was demonstrated that the carrier density can be 
controlled at will by electric gates. Therefore, if one can 
shrink the width of the single-layer graphite, the signa- 
ture of the novel ferromagnetism at nanoscale will start 
to emerge. 

In summary, combining the weak-coupling analysis and 
the DMRG technique, we have shown that the armchair 
NGR, due to electronic correlations and the flat band, 
has intrinsic magnetic moments. Furthermore, if gapless 
itinerant carriers are present, they mediate exchange cou- 
pling among the moments and lead to the ferromagnetic 
ground state. 
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